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IMPROVED  HANDLING-CHARACTERISTICS  MODEL 
I.  INTRODUCTION 

The  dynamics  of  motion  of  underwater  vehicles  is  basically  nonlinear  [l]-[4]. 
A major  thrust  of  this  year's  research  supported  by  NCSC  has  been  the  develop- 
ment of  nonlinear  models  of  handling  characteristics.  This  was  accomplished 
by  separating  the  steady-state  gain,  which  is  nonlinear,  and  the  dynamic 
transfer  function,  which  is  linear.  The  two,  placed  in  tandem,  model  the 
vehicle  behavior  over  wide  ranges  of  speed  and  control  surface  deflections. 
Further,  the  poles  of  this  transfer  function  more  accurately  reflect  the 
vehicle's  performance  than  those  of  the  traditionally  used  strictly  linear 
transfer  function.  For  example,  the  poles  of  the  traditional  transfer  function 
tend  to  show  greater  stability  than  the  vehicle  actually  possesses.  The  new 
model  overcomes  this  defect.  Indeed,  this  refinement  in  handling  qualities 
description  will,  it  is  hoped,  yield  a truer  correlation  between  pilot  rating 
of  the  craft,  the  handling  qualities  model  and,  when  desired,  the  simulator 
response  of  the  craft  for  pilot  training. 

The  nonlinear  steady-state  gain  is  obtained  by  determining  the  responses 
of  the  vehicle  to  various  degrees  of  step  inputs,  say  the  rudder  input,  and 
computing  the  ratio  of  the  final  response  value  to  the  input  magnitude.  The 
linear  dynamic  part,  i.e.,  the  transfer  function,  is  determined  by  the  GRAM 
Identifier,  a system  identification  technique  developed  in  the  ‘research  efforts 
of  1974-1976  [5]. 

We  will  develop  the  methodology  for  modeling  and  demonstrate  its  appli- 
cation to  a particular  submerged  vehicle  (USFRPV)  by  generating  the  vehicle 
responses  on  the  NCSC  trajectory  simulation  program  [5].  It  should  be  obvious, 
however,  that  the  ultimate  use  of  the  modeling  methodology  will  be  upon  tow 
tank  tests  of  actual  vehicles. 


II.  NONLINEAR  HANDLING-CHARACTERISTICS  MODELS 

To  find  a more  accurate  description  of  the  handling  characteristics  of  a 
vehicle  than  has  been  achieved  in  the  past,  we  propose  to  use  the  structure 
shown  in  Fig.  1.  Here  CNL  is  a static  nonlinearity,  custom  designed  to  fit 
the  vehicle's  actual  saturation  characteristic,  and  H(s)  is  a linear  transfer 
function  describing  the  transient  modes  of  the  vehicle.  Since  the  steady-state 
gain  of  H(s)  (namely  H(0))  is  a constant,  i.e.,  it  is  independent  of  the  ampli- 
tude of  the  input  to  H(s),  the  CNL  function  can  be  found  by  performing  a set 
of  step  function  tests  on  the  vehicle.  While  these  tests  can  of  course  be 
performed  in  a tow-tank  on  an  actual  vehicle,  we  will,  in  this  report,  use 
the  NCSC  trajectory  simulation  program  upon  a particular  test  vehicle,  the 
USFRPV.  Once  the  CNL  function  has  been  determined,  say  in  the  form  of  a table, 
it  is  shown  that  the  linear  transfer  function  H(s),  or  its  z-domain  equivalent 
H(z),  can  be  found  by  use  of  a linear  system  identification  program  developed 
under  a previous  research  contract  [ 5 ] , [ 6 ] . 


Fig.  1.  Nonlinear  model 


As  mentioned  above,  the  NCSC  Trajectory  Simulation  Program  is  utilized  here 
for  determining  the  responses  of  USFRPV.  The  following  changes  were  made  in 
the  program:  (i)  the  uu  terms,  such  as  Kuu  etc.,  were  zeroed  out,  and  (ii)  the 
vehicle  speed  was  synthetically  set  to  a constant  (by  setting  XD(1)  = 0). 

After  the  above  modifications  were  completed,  computer  runs  were  made  at 
speeds  of  5,  7.5,  10,  15,  20,  25,  30  knots  with  rudder  inputs  of  1,  5,  10,  20, 

35  degrees  at  each  speed.  These  runs  were  made  for  a 40  second  duration  with 
200  data  points.  The  rudder  input  shown  in  Fig.  2,  and  the  corresponding 
vehicle  responses  V,  W,  P,  Q,  R were  obtained  (and  punched  on  card  decks  for 
modeling  and  analysing) . 
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To  form  a CNL  table,  the  steady  state  value  of  response  Is  plotted  vs. 
the  corresponding  amplitude  of  step  Input.  Thus  a CNL  table  Is  obtained  for 
each  particular  speed  of  Interest.  This  process  Is  repeated  for  each  response 
variable.  Fig.  2 , for  example,  displays  the  steady  state  values  of  v for 
five  different  values  of  rudder  steps  (1°,  5°,  10°,  20°,  and  35°),  all  for  a 
vehicle  speed  of  5 knots.  A smooth  curve  through  these  points  yields  the 
CNL  characteristic  of  the  particular  response  variable  at  the  specified  speed. 


CNL 

Output 


Speed  (Knots) 

Fig.  2.  CNL  characteristics  curve  for  V versus  6^,  5 knots 

From  the  graph,  the  CNL  Output  (see  Fig.  2 ) corresponding  to  different  input 
amplitudes  over  a range  of  0 to  36°  can  be  read  and  tabulated.  CNL  tables  are 
prepared  in  this  manner  for  each  response  variable,  one  each  speed  of  interest. 

The  transfer  function  for  each  motion  variable  is  determined  via  GRAM 
identifier,  mentioned  previously,  using  the  CNL  table  and  the  input  and  output 
response  data.  The  input  pertains  to  the  rudder  whereas  the  output  corresponds  to 
V,  W,  P,  Q and  R.  The  CNL  value,  taken  from  the  CNL  table,  corresponding  to 
the  input  rudder  deflection  is  the  actual  input  to  the  transfer  function  block 
The  transfer  functions  used  for  the  response  variables  V,  W,  Q,  R are  of  the 
second  order.  However,  the  variable  P requires  a fourth  order  model  because 
of  visible  presense  of  a double  oscillation.  Transfer  function  models  were 
formed  for  all  variables  of  speeds  of  .7.5,  15,  25  knots;  a 5°  step  rudder  input 
as  used  in  each  run. 
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To  determine  the  effectiveness  of  the  Identified  model  a validation  phase 
was  included  in  our  investigation.  A Fortran  program,  VERIFY,  (see  Appendix  A) 
was  written  for  this  purpose.  It  predicts  the  vehicle  response  to  test  inputs 
and  compares  it  with  actual  (NCSC  TSP)  response.  Plots  exhibiting  this  com- 
parison as  well  as  a RMS  error  value  are  obtained.  Validation  runs  were  made 
at  (a)  the  speed  the  transfer  function  was  developed,  call  it  UQ,  (b)  a speed 
higher  than  Uq,  call  it  U+  and  (c)  a speed  lower  than  Uq,  call  it  U-. 

The  Z total  energy  between  the  vehicle  response  and  the  transfer  function 
response  are  in  Tables  1-3.  In  the  pages  following  the  tables,  the  CNL  curves, 
CNL  Tables,  and  the  validation  plots  are  given  for  each  speed  for  all  speed 
ranges  and  all  variables. 
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X OF  TOTAL  ENERGY 
5-10  SPEED  RANGE 


DEFLECTION 

— 

VARIABLE-SPEED 

1° 

5° 

O 

o 

20° 

35° 

39.87 

17.61 

4.82 

-6.47 

11.37 

V - 7.5  Kn 

26.83 

.144 

10.04 

17.17 

19.94 

16.29 

9.65 

17.03 

22.39 

24.27 

W - 5 Kn 

78.73 

37.94 

22.34 

5.40 

8.80 

W - 7.5  Kn 

38.63 

2.20 

17.32 

30.47 

37.01 

W - 10  Kn 

19.94 

23.48 

36.85 

33.99 

48.39 

P - 5 Kn 

111.0  . 

63.41 

24.81 

45.52 

54.84 

P - 7.5  Kn 

105.3  • 

6.88 

39.86 

60.48 

58.52 

P - 10  Kn 

98.85 

25.52 

49.12 

61.64 

47.88 

Q - 5 Kn 

111.9 

57.67 

19.14 

22.62 

31.44 

Q - 7.5  Kn 

83.61 

3.99 

31.31 

43.50 

40.37 

Q - 10  Kn 

55.05 

24.22 

40.60 

45.11 

42.08 

R - 5 Kn 

35.56 

12.89 

1.42 

8.89 

12.77 

R - 7.5  Kn 

21.02 

0.47 

9.33 

15.81 

18.17 

R - 10  Kn 

17.52 

6.99 

13.99 

19.16 

20.94 

Table  2 


2 OF  TOTAL  ENERGY 
10-20  SPEED  RANGE 


VARIABLE-SPEED 


- 10  Kn 

- 15  Kn 

- 20  Kn 


W - 10  Kn 
W - 15  Kn 
W - 20  Kn 


- 10  Kn 

- 15  Kn 

- 20  Kn 


DEFLECTION 

1° 

5° 

10° 

20° 

35° 

37.94 

11.05 

3.26 

3.73 

6.25 

21.62 

.412 

6.07 

11.07 

13.83 

11.28 

6.46 

11.54 

15.71 

21.15 

94.99 

35.17 

15.80 

5.38 

6.97 

47.33 

2.44 

11.63 

19.20 

21.36 

19.19 

14.45 

21.99 

26.67 

26.26 

137.5 

23.33 

32.50 

42.27 

24.59 

95.84 

3.30 

17.56 

15.76 

165.92 

84.92 

31.19 

49.85 

161.0 

137.5 

104.8 

21.36 

15.63 

11.31 

14.71 

59.34 

3.52 

14.45 

27.19 

42.35 

41.63 

20.07 

30.68 

40.70 

50.20 

33.68 

8.48 

1.51 

5.24 

7.54 

17.41 

.716 

5.57 

9.97 

12.06 

12.63 

4.94 

9.31 

13.34 

15.92 

Table  3 

Z OF  TOTAL  ENERGY 
20-30  SPEED  RANGE 


DEFLECTION 


VARIABLE-SPEED 


- 20  Kn  21.89 

- 25  Kn  15.36 

- 30  Kn  10.69 


W - 20  Kn  AO. 29 
W - 25  Kn  24.82 
W - 30  Kn  17.94 


P - 20  Kn  43.44 
P - 25  Kn  35.60 
P - 30  Kn  36.12 


Q - 20  Kn  33.07 
Q - 25  Kn  23.70 
Q - 30  Kn  27.66 


R - 20  Kn  19.42 
R - 25  Kn  13.86 
R - 30  Kn  10.44 


5 


4.99 

.405 

3.43 


37.79 

2.86 

59.93 


11.30 

2.93 

7.73 


3.43 

.798 

2.67 


o 

O 

iH 

20° 

35° 

1.24 

5.17 

11.60 

4.88 

10.11 

16.70 

8.17 

15.35 

17.96 

2.09 

6.84 

13.12 

7.65 

13.88 

13.60 

12.94 

15.70 

19.50 

28.74 

66.59 

110.8 

62.35 

113.9 

107.5 

111.4 

110.0 

106.5 

4.92 

8.32 

18.31 

8.50 

15.50 

13.61 

14.65 

15.41 

18.77 

1.80 

6.25 

9.36 

4.66 

9.25 

10.49 

7.58 

11.54 



11.32 

Z - 1.950Z  + 9.512  x 10 


Poles  in  the  s Domain* 

-1.25  x 10_1  + j 8.49  X 10-2 


♦For  conversion  from  Z-domain  Co  S-domaln,  see  references  [11], [8], 


7.5,  10  knots  V/6r  CNL  curves 


Input  Deflection  Magnitude  (In  degrees) 
Fig.  3 
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15,  20,  25,  and  30  knot  V/6^  CNL  curves 


Input  Deflection  Magnitude  (In  degrees) 


Table  5 

Custom  nonlinearity  (CNL) : V versus 
10-20  knot  range 


s 

r 

V steady  state 

Degree 

10  knots 

15  knots 

20  knots 

0 

0.000 

0.00 

0.00 

1 

0.143 

0.22 

0.29 

2 

0.215 

0.32 

0.43 

3 

0.270 

0.40 

0.54 

4 

0.316 

0.47 

0.62 

5 

0.357 

0.53 

0.70 

6 

0.391 

0.59 

0.77 

7 

0.427 

0.64 

0.83 

8 

0.455 

0.68 

0.89 

9 

0.483 

0.72 

0.94 

10 

0.509 

0.76 

0.99 

11 

0.533 

0.80 

1.04 

12 

0.555 

0.83 

1.08 

13 

0.579 

0.86 

1.12 

14 

0.600 

0.89 

1.16 

15 

0.620 

0.92 

1.20 

16 

0.640 

0.95 

1.23 

17 

0.660 

0.98 

1.26 

18 

0.680 

1.00 

1.29 

19 

0.700 

1.03 

1.33 

20 

0.720 

1.06 

1.36 

21 

0.739 

1.09 

1.38 

22 

0.755 

1.12 

1.40 

23 

0.771 

1.14 

1.42 

24 

0.790 

1.16 

1.44 

25 

0.809 

1.18 

1.46 

26 

0.823 

1.21 

1.47 

27 

0.841 

1.23 

1.48 

28 

0.860 

1.25 

1.50 

29 

0.877 

1.27 

1.52 

30 

0.893 

1.29 

1.53 

31 

0.912 

1.31 

1.54 

32 

0.930 

1.33 

1.56 

33 

0.946 

1.34 

1.57 

34 

0.962 

1.36 

1.58 

35 

0.980 

1.38 

1.59 

36 

0.999 

1.40 

1.60 

Fig.  8.  V Versus  6^,  10  Knots 
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Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  v/sr 

Model: 


25  Knots 


U_!  20  Knots 

U+:  3Q  Knots 

H6r(Z)  Determined  Using  U0  and  5°  «r  Step 

r 7 x 10'2z  - 2.615  x 10'° 


Z - 1.854Z  + 8,630 

Poles  in  the  s Domain* 

-3,683  x 10"1  + j 3,439  x 10_1 


*For  conversion  from  Z-domain  to  S— domain,  see  references  [11], [8] 
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Table  6 

Custom  nonlinearity  (CNL) : V versus  6 ^ 
20-30  knot  range 

V steady  state 


SiSiSiSSSti 


Results  of  Validation  Tests  For: 

USF  RPV 

Transfer  Function:  w/6r 

1 3 

Model: 

* . • 1.  i 

• • I 

i | 

TMI 

r 

ll5r(Z)  

V 

• "_,im 

7.5  Knots 

; 

U.: 

5 Knot's 

1 

. 

U+! 

10  Knots 

. I 

‘ 

ij 

H*r(Z)  Determined  Using  Uq  and  5 

# «r  Step 

' 

Hfr  * . 

-2.562  x 10~3Z  + 2.707 

■ 

-3 

x 10 

-2.002Z  + 1.003 

Poles  in  the  s Domain* 

-3  * -1 

7.071  x 10  J + J 9.211  x 10 

•For  conversion  from  Z-domain  to  S-domain,  see  references  [11], [8] 
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ONLn^WNJHOvOC»NjCJ'Ln^UN)HOvOODv4(J\Ln^UJNHOv003vjO\Oi^UK)HO 


<5 

r 

Degree 


Table  7 

Custom  nonlinearity  (CNL) : w versus  6^ 


5-10  knot  range 

w steady  state 


5 knots 


7.5  knots 


10  knots 


0.000 

0.000 

0.000 

0.004 

0.013 

0.026 

0.015 

0.035 

0.057 

0.023 

0.054 

0.060 

0.035 

0.073 

0.112 

0.043 

0.096 

0.135 

0.057 

0.113 

0.156 

0.067 

0.127 

0.175 

0.076 

0.142 

0.190 

0.086 

0.155 

0.207 

0.094 

0.166 

0.219 

0.105 

0.176 

0.235 

0.115 

0.183 

0.249 

0.122 

0.195 

0.262 

0.130 

0.206 

0.275 

0.137 

0.215 

0.290 

0.145 

0.223 

0.302 

0.150 

0.231 

0.313 

0.157 

0.240 

0.323 

0.162 

0.247 

0.333 

0.165 

0.257 

0.340 

0.170 

0.264 

0.349 

0.172 

0.267 

0.355 

0.175 

0.273 

0.361 

0.177 

0.277 

0.367 

0.180 

0.281 

0.373 

0.183 

0.285 

0.377 

0.185 

0.288 

0.383 

0.187 

. 0.290 

0.387 

0.189 

0.294 

0.393 

0.191 

0.296 

0.398 

0.193 

0.300 

0.402 

0.195 

0.303 

0.406 

0.197 

0.305 

0.412 

0.198 

0.306 

0.415 

0.203 

0.307 

0.420 

0.203 

0.310 

0.423 

27 


txi 


ftt  5 TO  1*  KKST  1TOS0F  V/ BEL  9 


10  Ot&fEf  S KKOT5 


Fig.  15.  W Versus  5 Knots 
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Fig.  16.  W Versus  6^,  7.5  Knots 
29 


Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  w/6r 

Model: 


UqS  15  Knots 

U_!  10  Knots 

U+:  20  Knots 

H?,r(Z)  Determined  Using  U.  and  5°  sr  Step 

W u 

H^r  * -1.046  x 10~2Z  + 1.146  x 10~2 

Z2  -1.966Z  + 9.667  x 10_1 

Poles  in  the  s Domain* 

-8i.  460  x 10~2  + j 1.356  x 10_1 


*For  conversion  from  Z-domain  to  S-domain,  see  references  [11] >[8] 
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15,  20,  25,  30  knots  W/6r  CNL  curves 


lef  lect 


iifiiifiiiil 


10  DEORFF.  15  KNOTS 

Fig.  20.  W Versus  6^,  15  Knots 
35 


Itf  14  111  Ml  Ml  Ml  M 

10  OF GREF  20  KNOTS 


Fig.  21.  W Versus  6^,  20  Knots 


1: 
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Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  w/6r 

Model: 


U 


o* 


25  Knots 


H: 

U 


20  Knots 


30  Knots 


H6r (Z)  Determined  Using  Un  and  5°  Step 


J|5r  , -1.191  x 10~2Z  + 1.467  x 10~2 

W ZZ  -1.931Z  + 9.337  x 10_1 


Poles  in  the  s Domain' 


-1.716  x 10_1  + j 1.992  X 10* 1 


•For  conversion  from  Z-domain  to  S-domain,  see  references  [11] ,[8]. 
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TABLE  9 

Custom  Nonlinearity  (CNL) : W Versus  6r 
20-35  Knot  Range 


r 

agree 


R Steady  State 
20  knots  25  knots 


30  knots 


0.00 

0.07 

0.13 

0.19 

0.23 

0.28 

0.32 

0.37 

0.40 

0.44 

0.48 

0.52 

0.55 

0.58 

0.61 

0.65 

0.68 

0.71 

0.74 

0.78 

0.81 

0.85 

0.88 

0.92 

0.96 

0.99 

1.03 

1.06 

1.10 

1.14 

1.18 

1.21 

1.24 

1.28 

1.32 

1.35 

1.38 


f 


st  n m oars  *mcf  « / oti  t 


20  to  so  dot  wwof  m r oei  » 


20  to  38  KNcrs  « t oei  r 


Fig.  23,  W Versus  5^,  25  Knots 
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Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  p/6r 

Model: 


Uo*  7.5  Knots 

(J  ; 5 Knots 

(J^j  10  Knots 

H6r(Z)  Determined  Using  Ud  and  5°  Step 

p 

fj5r  - -2.633  x lO^Z3  + 8.342  x lO^Z2  - 8.796  x lO^Z  + 3.086  x 10 
Z4  -3.837Z3  + 5.589Z2  - 3.669Z  + 9.165  x 10_1 

Poles  in  the  s Domain* 

4.3  x 10"2  + j 5.411  X 10~2 

*For  conversion  from  Z-domain  to  S-domain,  see  references  [11], [8]. 
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Table  10 


6 

r 

Degree 


Custom  nonlinearity  (CNL) : P versus  6 

r 

5-10  knot  range 

| P steady  state 


5 knots 

7.5  knots 

0.00 

0.0 

-0.01 

-0.03 

-0.02 

-0.06 

-0.03 

-0.09 

-0.045 

-0.12 

-0.06 

-0.16 

-0.072 

-0.20 

-0.084 

-0.23 

-0.096 

-0.25 

-0.108 

-0.28 

-0.120 

-0.30 

-0.129 

-0.33 

-0.138 

-0.36 

-0.147 

-0.38 

-0.156 

-0.40 

-0.165 

-0.42 

-0.174 

-0.44 

-0.183 

-0.46 

-0.193 

-0.48 

-0.201 

-0.50 

-0.210 

-0.52 

-0.217 

-0.54 

-0.224 

-0.56 

-0.231 

-0.57 

-0.236 

-0.58 

-0.240 

-0.60 

-0.243 

-0.61 

-0.246 

-0.62 

-0.249 

-0.63 

-0.252 

-0.64 

-0.255 

-0.65 

-0.258 

-0.66 

-0.261 

-0.67 

-0.264 

-0.68 

-0.267 

-0.70 

-0.270 

-0.71 

-0.273 

10  knots 

0.0 

-0.07 

-0.13 

-0.19 

-0.25 

-0.31 

-0.37 

-0.43 

-0.49 

-0.54 

-0.59 

-0.64 

-0.68 

-0.72 

-0.76 

-0.81 

-0.86 

-0.90 

-0.95 

-0.99 

-1.04 

-1.07 

-1.10 

-1.13 

-1.16 

-1.19 

-1.22 

-1.25 

-1.28 

-1.31 

-1.33 

-1.35 

-1.37 

-1.39 

-1.41 

-1.42 

-1.44 
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Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  p/$r 

Model: 


Uc : 15  Knots 

U_ ; 10  Knots 


20  Knots 

Hpr(Z)  Determined  Using  Ud  and  5°  Step 

H6r  - -1-604  x lo'V  + 6.043  x IO'1!2  - 6.702  x IQ*1  + 2.464  x IQ-1 
Z4  - 3.771Z3  + 5.370Z2  - 3.430Z  + 8.305  x 10_1 

Poles  in  the  s Domain* 

1-99  x io"1 
. -5.9985  x 1.0'2 
-5.34  x IO-1  + j 1.244 

*For  conversion  from  Z-domaln  to  S-domaln,  see  references  [11] , [8]. 
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Table  11 

Custom  nonlinearity  (CNL) : P versus  6^ 
10-20  knot  range 


Fig*  29.  P Versus  6^,  15  Knots 


j 
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Results  of  Validation  Tests  For:  USF  RPV 


Transfer  Function: 


Model: 


CNLi 


H*r(Z) 


25  Knots 


20  Knots 


30  Knots 


Hp(Z)  Determined  Using  Uq  and  58  «r  Step 


l6r  - -1.240  x lO^Z3  + 4.374  x lO'V  - 5.026  x IQ-1 

Z4  -3.641Z3  + 4.993Z2  - 3.060Z  + 7.071  x 10_1 


Poles  in  the  s Domain* 


-4.699  X 10 


-3.665  x 10 


♦For  conversion  from  Z-domain  to  S-domain,  see  references  [11], [8], 
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Table- 12 


its  * 


n?m 


T 

* 


Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  q/6t 

Model: 


U 


o* 


U_: 

U+! 


7.5  Knots 


5 Knots 


10  Knots 


H6  (Z)  Determined  Using  U.  and  5°  «r  Step 
q 


H6r  - -9.526  x 10~2  Z + 1.006  x 10~2 

q z2  - 1.994Z  + 9.945  X 10_1 

Poles  in  the  s Domain* 


-1.375  x 10-2  + j 1.258  x 10_1 


*For  conversion  from  Z-domain  to  S-domain,  see  references  [11] ,[8], 
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7.5,  10  knots  Q/6  CNL  curves 


Input  Deflection  Magnitude  (In  degrees) 
Fig.  33 


a>U1>WKJHO'O00vJO'W*.UMHO'O06vJ»UlMJIOHO'O05'Ja'Ult'WS)HO 
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Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  Q/5r 

Model: 


II  • 15  Knots 

*'o* 

U_!  10  Knots 

U+:  20  Knots 

H*r(Z)  Determined  Using  Ud  and  5°  $r  Step 

H6r  - -1,281  x 10~2Z  + 1.453  x 1 0~2 

q Z2  - 1.974Z  + 9.754  x 10_1 

Poles  in  the  s Domain* 

-6.219  x 10-2  + j 1.758  x 10_1 


*For  conversion  from  Z-domain  to  S-domain,  see  references  [11], [8]. 
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1 


15,  20,  25,  30  knots  £/6  CNL  curves 


Deflection  Magnitude  (In  degrees) 
Fig.  38 


10  OEGRtF  10  KNOTS 


Fig.  39.  Q Versus  5^,  10  Knots 
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. -jrzrrxz 


10  TO  » KOTS  KMCf  0 I *i  t 


M «4l  111  Ml  It  i * t 

35  DEGREE  15  KNOTS 


10  DEGREE  15  KNOTS 


Fig.  40.  Q Versus  6r>  15  Knots 
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iai  lit  li  aV  AY  A.  AY  AY 

1 DECREE  20  KNOTS 

a • 

» t 

i 

Mt 

ta 

I.N 

»»  10  TO  20  KNOTS  KflNOf  C / OEL  » 

L * 1 * * »-» 

, 

4. 

Ml 

i.a 

* *•  LJ  111  |«  t M|  Ml 

™ nFGREE  20  KNOTS 


10  ro  ?0  KNOTS  «*G€  Q / DO.  f 


tii  ill  t«4  at  i mi  aii  at  it  i n: 

5 DECREE  20  KNOTS 


10  TO  20  KNOTS  *WG£  O / OEL  0 


. 1 .... L ....  1 


*•  l*  Ml  HI  It  1 Ml 

95  DECREE  20  KNOTS 


10  DECREE  20  KNOTS 


Pig.  41,  0 Versus  5r,  20  Knots 
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Results  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  q/&t 

Model: 


Uo:  25  Knots 

U_ : 20  Knots 


u. 


30  Knots 


H?r(Z)  Determined  Using  U„  and  5*  *r  Step 


H6r  - -5.831  x 10~2Z  + 8.861  x 10~2 

q Z2  - 1.939Z  + 9.423  x 10"1 

Poles  in  the  s Domain* 

-1.486  x 10"1  + j 2.283  x 10_1 


*For  conversion  from  Z— domain  to  S— domain,  see  references  [11], [8], 
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m \ m 14  hi  i% . «• » * • » . ii  » « «f 

10  KCPEF  30  •'MOTS 


Fig.  44.  Q Versus  30  Knots 


7.5,  10  knot  R/6  CNL  curve 


Table  16 


Custom  nonlinearity  (CNL) : R versus  6r 
5-10  knot  range 


6 

r 

R steady  state 

Degree 

5 knots 

7.5  knots 

10  knots 

0 

0.0 

0.00 

0.00 

1 

-0.32 

-0.53 

-0.79 

2 

- .55 

- .90 

-1.20 

3 

- .72 

-1.15 

-1.55 

4 

- .88 

-1.38 

-1.85 

5 

-1.06 

-1.59 

-2.10 

6 

-1.16 

-1.77 

-2.35 

7 

-1.30 

-1.94 

-2.56 

8 

-1.40 

-2.08 

-2.77 

9 

-1.50 

-2.22 

-2.95 

10 

-1.60 

-2.35 

-3.13 

11 

-1.67 

-2.47 

-3.30 

12 

-1.77 

-2.60 

-3.47 

13 

-1.85 

-2.71 

-3.63 

14 

-1*93 

-2.83 

-3.78 

15 

-2.00 

-2.95 

-3.95 

16 

-2.07 

-3.05 

-4.08 

17 

-2.15 

-3.17 

-4.23 

18 

-2.23 

-3.27 

-4.38 

19 

-2.30 

-3.40 

-4.52 

20 

-2.36 

-3.52 

-4.68 

21 

-2.44 

-3.60 

-4.80 

22 

-2.50 

-3.70 

-4.95 

23 

-2.57 

-3.80 

-5.08 

24 

-2.63 

-3.92 

-5.23 

25 

-2.70 

-4.02 

-5.35 

26 

-2.78 

-4.13 

-5.55 

27 

-2.84 

-4.23 

-5.63 

28 

-2.90 

-4.32 

-5.77 

29 

-2.96 

-4.44 

-5.90 

30 

-3.04 

-4.55 

-6.04 

31 

-3.10 

-4.65 

-6.16 

32 

-3.16 

-4.75 

-6.30 

33 

-3.23 

-4.85 

-6.44 

34 

-3.33 

-5.00 

-6.58 

35 

-3.37 

-5.06 

-6.70 

36 

-3.43 

-5. 

-6.85 
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10  KGRCf  5 r.MiTS 


Pig.  46.  R Versus  5 Knots 
77 


♦ ti  U in  •«  I n l HI  Ml  B 1 »|  HI 

35  DECREE  7.5  KNOTS 


l m u ill  itt  r i hi  a. 

JO  DECREE  7.5  KNOTS 


Fig.  47.  R Versus  6r>  7.5  Knots 
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??? 


10  DtGfiEF  10  KNOTS 


Fig.  48.  R Versus  6r>  10  Knots 
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80 


15,  20,  25,  and  30  Knot  R/6r  CNL  curves 


Deflection  Magnitude  (In  degrees) 
Fig.  49 


(f  iisiiii 


Fig.  52.  R Versus  6^,  20  Knots 
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JIesults  of  Validation  Tests  For:  USF  RPV 

Transfer  Function:  R/6r 

Model: 


Uo:  25  Knots 

U_:  20  Knots 

U+:  30  Knots 

H*r(Z)  Determined  Using  Ud  and  5“  «r  Step 

H6r  - 3.557  X 10~2Z  - 2.615  x 10~2 

r z2  -1.854Z  + 8.630  x 10_1 

Poles  in  the  s Domain* 

-4.223  X 10" 1 + j 2.980  x 10_1 


*Tot  conversion  from  Z-domain  to  S— domain,  see  references  [11], [8], 
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10  CEUREf  30  KNOTS 


Fig,  55,  R Versus  6^,  30  Knots 
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III.  AMPLITUDE  AND  SPEED  NORMALIZATION 

In  the  equations  of  motion,  the  force  term  produced  by  a rudder  deflection 

6 is  proportionate  to  6 * U where  U is  the  vehicle  speed.  It  is  therefore 

r r o o 

useful  to  normalize  vehicle  responses  by  a factor  C ■ 6 * U where  U is  in  ft. 

r o o 

per  second.  Such  normalized  responses  are  given  in  this  section  for  the  y and  V 
variables.  Fig.  56  to  58,  Fig,  62  to  65  present  a comparison  of  Y and  W responses, 
respectively,  to  various  amplitude  steps,  each  figure  drawn  for  a fixed  speed. 

Also,  it  is  known  from  principles  of  hydrodynamics  that  the  responses  of 
a submerged  vehicle  involves  a time-scale  factor  which  is  directly  proportional 
to  speed.  Specifically,  this  means  the  following.  If  the  transfer  function  of 
a particular  response  variable,  say  side-slip,  is  H(s^)  at  a forward  speed  Uq^» 
and  it  is  H^)  at  another  value  of  forward  speed,  » then 


Clearly,  a time  normalization  must  be  performed  if  the  transient  behavior  at  two 
different  speeds  are  to  be  compared.  This,  of  course,  will  be  in  addition  to 
the  amplitude  normalization  discussed  above. 

Fig . 59  to  61,  Fig.  66  to  68  give  comparison  of  amplitude  and  time 

normalized  V and  W,  respectively,  responses  at  various  speeds,  each  figure  drawn 
for  a fixed  amplitude  rudder  step. 


91 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
PROM  COPY  FURNISHED  TO  ICC  — 


1.1---., 


7- 5 KM 


DT~0 . 2 

— DT=0 . 2k 
— DT=0.30 
— DT=0.A0 
— DT=0.60 
— OT=O.SO 
DT=1 . 20 


■ i~™  rrr'r-rr^T^~^-Tr"  - ■ i 

too 


Key  To  Constant  Angle  Plots 


Fig.  59 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  EURKISHJSD  TO  DOC 


WY  KMSTRNT  SpEf:D  W /Cul.  R 


IV.  CNL  NORMALIZATION  OF  RESPONSES 


It  was  demonstrated  in  Section  II  that  the  vehicle  responses  bear  a nonlinear 
relationship  to  the  amplitude  of  the  rudder  input.  In  particular  the  steady  state 


gain,  for  example  Hv(0)  = for  step-slide  to  a step  rudder,  exhibits  a 


saturation  behavior.  In  order  to  compare  the  transient  modes  at  various  rudder 
amplitudes,  we  must  remove  the  effect  of  nonlinear  gain  H(0).  Thus,  we  must  deter- 
mine CNL  normalized  responses  by  dividing  the  amplitude  normalized  response  by  a 
CNL  gain  reduction  factor  at  the  particular  rudder  deflection.  For  example. 


r 


v(t) 


f (6  )*6  *Un 
r r 0 


where  v(t)  denotes  the  CNL  normalized  response  and 
CNL  gain  at  6 


^^r^  CNL  gain  at  1°  rudder 


Also,  of  course,  time  normalization  as  discussed  in  Section  III  (p.  91) 
must  be  performed  if  transient  response  at  different  speeds  are  to  be  compared. 


That  is,  if  v(t^)  and  v(t^)  represent  V and  W 


responses  at  two  different 


vehicle  speeds,  Uq1  and  U £>  respectively,  then 


fc2  = 


01 


02 


Fig.  69  to  72,  Fig.  73  to  76  illustrate  the  effects  of  CNL  Normalization 
on  the  transient  response  for  V and  W responses,  respectively, 
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V.  PARAMETRIZATION  OF  TRANSFER  FDNCTIONS 


The  undimensionalized  force/moment  input  in  the  equation  of  motion  is 

y - 5 * U 
r o 

so  that  the  vehicle  characterization  can  be  expected  to  be  different  for 
various  values  of  y (in  spite  of  the  use  of  CNL).  Even  for  a fixed  speed 
the  characterization,  in  particular  H(s),  can  be  expected  to  deviate  from 
one  nominal  function.  To  surmount  this  problem  we  propose  to  parameterize 
the  coefficients  of  H in  terms  of  y;  or,  at  a fixed  speed,  in  terms  of  6^. 
In  this  section  we  will  deal  with  constant  speed  characterization. 


To  advoid  this  problem  of  incompatible  transfer  functions,  a different 
strategy  was  needed  in  which  the  coefficients  of  the  basic  second  order  transfer 
functions  were  varied.  Advoiding  the  need  to  generate  transfer  functions  at 
all  rudder  inputs  applied  to  the  vehicle,  we  chose  transfer  function  coefficient 
parametrlzation  as  the  practical  approach  to  the  problem.  From  the  basic 
transfer  function  configuration 
cs  + b*  , 


H(s) 


8 + as  + b 


it  can  be  observed  that  we  need  to  parametrize  a,  b and  c in  te’frms  of  6^.  (From 
the  transfer  functions  in  Tables  and  , it  can  be  seen  that  b*  is 
approximately  b,  thereby  reducing  the  number  of  equations  required). 

The  parameterization  equations  were  determined  for  response  variables  V 
and  W at  speeds  of  7.5,  15  and  25  knots.  At  each  speed,  s domain  transfer 
functions  were  determined  with  step  rudder  inputs  of  1°,  5°,  10°,  20°  and  35°. 
These  transfer  functions  are  shown  in  Tables  19  and  20  .We  decided  to  parametrize 
a,  b and  c in  the  form 

£ » a + 0 (6  ) + y (5  )2 
a a r 'a  r 


“b + W + 


a + BP(«J  + vA6ry 

c c r c r 


The  coefficients,  a,  0,  y are  determined  using  the  least  squares  technique. 
Consider,  for  example,  the  parametrlzation  of  a.  The  equations  are: 
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d 


*1° 

a + 6(1°)  + y (1°)2 

a5° 

a + 3(5°)  + Y (5s)2 

£ 

O 

1 

a + 3(10°)  + Y(10°) 

a20°  " 

a + 3(20°)  + y(20°) 

a35°  " 

a + 3(35°)  + y(35°) 

or  In  matrix  form 
a - X e 

T 

where  0^  = [a,  3,  yl  • By  the  least  squares  technique, [10] 

0 - (XT  W-1  X)_1  X W-1  a 

where  W is  the  weighting  matrix  W = Diag{l,  4,  9,  16,  25} 

In  Table  21  thru  26 , the  parametrization  equations  are  given  along  with 
a comparison  of  parametrized  coefficients  (a,  b,  c)  to  the  transfer  function 
coefficients  at  various  rudder  deflections. 


i 
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V Transfer  Functions 


W Transfer  Functions 


Parameterization  Equations:  r 


a - 0.2091  + 0.00896r  - 0.000046r2 

b - 0.0057  + 0.0034Sr  - 0.000026r2 

c - 0.3275  - 0.03876r  + 0.000946r2 

Comparison  of  Coefficients: 

a - transfer  function  coefficient 


a - parameterized  coefficient 


<5r 

a 

A 

a 

b 

6 

c 

A 

c 

1° 

0.307 

0.290 

0.009 

0.009 

0.218 

0.218 

5° 

0.052 

0.158 

0.023 

0.023 

0.250 

0.253 

10° 

0.061 

0.034 

0.04 

0.039 

0.307 

0.295 

O 

O 

CM 

0.069 

-0.071 

0.065 

0.068 

0.365 

0.375 

35° 

0.057 

0.122 

0.107 

0.106 

0.491 

0.487 

1?8 


a - 0.27491  + 0.04141<5r  - 0.000946r2 
b - 0.00644  + 0.020386r  - 0.000396r2 
c - 0.06162  + 0.002625r  - O.OOOOlfir2 


Comparison  of  Coefficients: 


a - transfer  function  coefficient 
a - parameterized  coefficient 


6r 

a 

A 

a 

b 

b 

c 

c 

1° 

0.310 

0.315 

0.027 

0.0264 

0.06 

0.0642 

5° 

0.489 

0.458 

0.096 

0.0987 

0.098 

0.0745 

10° 

0.596 

0.595 

0.196 

0.071 

0.09 

0.0868 

20® 

0.673 

0.726 

0.267 

0.259 

0.066 

0.110 

35® 

0.593 

0.570 

0.242 

0.245 

0.160 

0.141 
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TABLE  23 

Variable : V 

Speed:  25  Knots 

Transfer  Function  in'  the  Form: 

H(.)  ■ ■ * b 

s + as  + b 

Parameterization  Equations 

a - 0.48255  + 0.036466r  - 0.001135r2 

f>  - 0.5395  + 0.03376r  - 0.000976r2 

c - 0.08347  + 0.0171l6r  - 0.000076r2 
Comparison  of  Coefficients 

a*-  transfer  function  coefficient 
a - parameterized  coefficient 


6r 

a 

A 

a 

b 

A 

b 

c 

/V 

c 

1° 

0.500 

0.518 

0.0750 

0.0867 

0.103 

0.101 

5° 

0.737 

0.637 

0.254 

0.198 

0.164 

0.167 

O 

O 

rH 

0.743 

0.734 

0.340 

0.294 

0.204 

0.248 

• 

O 

eg 

0.577 

0.759 

0.182 

0.340 

0.467 

0.399 

35° 

0.450 

0.372 

0.110 

0.0465 

0.578 

0.602 

TABLE  24 

Variable : W 

Speed : 7.5  knots 

Transfer  Function  in  the  Form 

H(s)  - ■+-  > 

s + as  + b 

Parameterization  Equations 


a - -0.11588  + 0.034156r  - 0.000816* 


b - 0.00721  + 0.00033<5r  + 06 

r 

# 

c - -0.0028  - 0.002336r  + 0.000036" 

l 

Comparison  of  Coefficients 

a - transfer  function  coefficient 
a - parameterized  coefficient 


6 

r 

a 

a 

b 

8 

c 

A 

c 

1° 

-0.0814 

-0.0825 

0.0076 

0.00754 

-0.00458 

-0.00458 

5° 

-0.0141 

0.0346 

0.0085 

0.00887 

-0.0132 

-0.0132 

10° 

0.315 

0.144 

0.0106 

0.0105 

-0.023 

-0.0228 

20° 

0.0851 

0.242 

0.0145 

0.0140 

-0.0374 

-0.0376 

35° 

0.133 

0.0852 

0.019 

0.0192 

-0.0494 

-0.0493 

TABLE  25 


Variable:  W 

Speed : 15  knots 

Transfer  Function  in  the  Form: 

H(s)  - -&-+■ b 

s + as  + b 

Parameterization  Equations: 

a - 0.02412  + 0.023236r  - 0.000326r 

b - -0.0031  + 0.0056r  - 0.000162 

r 

c - -0.0102  - 0.0086r  + 0.000252 

r 

» 

Comparison  of  Coefficients: 

a - transfer  function  coefficient 
a - parameterized  coefficient 


a 

A 

a 

b 

A 

b 

c 

A 

c 

0.0405 

0.0470 

0.000904 

0.0018 

-0.0163 

-0.018 

0.169 

0.1322 

0.0255 

0.0202 

-0.0557 

-0.0463 

0.227 

0.2242 

0.0399 

0.0398 

-0.076 

-0.0745 

0.294 

0.3599 

0.0593 

0.0683 

-0.0896 

-0.108 

0.471 

0.4426 

0.0883 

0.0843 

-0.1062 

-0.0985 

TABLE  26 


Variable:  W 

Speed:  25  Knots 

Transfer  Function  in  the  Form: 

„ , s cs  + b 

H(s)  = —= 

s + as  + b 

Parameterization  Equations: 

a = 0.1300  + 0.03535r  + 0.001352 

b = 0.00648  + 0.000626r  + 0.0006462 

r 

c = -0.03133  - 0.006046r  -0.0003952 

r 

Comparison  of  Coefficients: 

a - transfer  function  coefficient 


A 

a - 

parameterized 

coefficient 

6r 

a 

a 

b 

A 

b 

c 

c 

1° 

0.169 

0.167 

0.0222 

0.00775 

-0.0394 

-0.0378 

5° 

0.343 

0.338 

-0.0691 

0.0257 

-0.0685 

-0.0713 

h-4 

O 

o 

0.534 

0.610 

0.125 

0.0771 

-0.105 

-0.131 

20° 

1.44 

1.34 

0.369 

0.277 

-0.349 

-0.308 

35° 

2.88 

2.91 

0.771 

0.817 

-0.706 

-0.721 
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00010 

00020 

00030 

00040 

00050 

00060 

00070 

00080 

00090 

001O0 

00110 

00120 

00130 

00140 

00150 

00160 

00170 

00180 

00190 

00200 

00210 

00220 

00230 

00240 

00250 

00260 

00270 

00280 

00290 

00300 

00310 

00320 

00330 

00340 

00350 

00360 

00370 

00380 

00390 

00400 

00410 

00420 

00430 

00440 

00450 

00460 

00470 

00480 

00490 

00500 

00S10 

00520 

00530 

00540 

00550 


//STEPZ  EXEC  FTG1CLG,LIB5-'USF.  IR.PLOTLIB' 

//FORT . SYS  I N DD  * 

C • 

C 

C PROGRAM  VERIFY 

C SIMULATES  A SINGLE  INPUT  - SINGLE  OUTPUT  CNL  ♦ H(Z>  CASCODE 
C COMPARES  THE  RESPONSE  WITH  NCSL  PROGRAM  TRAJECTORY 


DIMENSION  X(200,15),Y(200,15),YREC(200,1S),YPLOT(200,3), 
A NAMEM(1S,15),NMAI0(15,15),UNL(  37 ) ,BUFF ( 200, 4 ) , SP< 10 ) , 

B TITLEC80) 

REAL»8  GAMMA (20) ,DELTA,GAI N,GAMMA1(20, 8 ) 

REAL *8  XLAMDA( 20) 

WR I TE( 6,202 ) 

READ(5,l)N,NMOD, I XX, I YY, DELTA 
WRITE(6,11) NMOD, N 
WRITE(G,203) 

NP1-N+1 

NP2«N*2 

NPNP1-N+N+1 

NPNP2*N+N*2 


READ  ALL  H(Z)  MODELS  ( E. 


FOR  THREE  DIFFERENT  SPEEDS,  10,15,20 


DO  199  J "l, NMOD 

READ ( 5 , 015 ) ( NAMEM( I ,J),  1-1,15) 

READ(5,2)  ( GAMMA (I ) , 1 «1 , NP1 ) 

HRI TE ( 6, 4 ) 

CALL  PRV EC ( GAMMA ( 1 ) ,NP1) 

READ(  5, 2 ) ( GAMMA ( I ) , I «=NP2,NPNP2) 

WRITE(6,5) 

CALL  PRVEC(GAMMA(NP2 ) , NP1 ) 

D0546I*NP2,NPNP2 
GAMMA ( I ) —GAMMA (I  ) 

DO  199  I »1,NPMP2 
GAMMA1  ( I , J ) *=  GAMMA  ( I ) 

READ  CNL  AND  PLOT  OPTIONS,  AND  CNL  TABLE 
READ ( 5 , 1 ) I CNL  ,IPLOT,IDYN 
WR I TE ( 6, 12 ) I CNL, I PLOT, I DYN 
I F ( I CNL. NE.0)READ(5, 9) (UNL(K) , K»1 ,37) 

READ  NUMBER  OF  MANUEVERS  (FOR  EACH  MODEL), THEIR  TITLES, 

THEIR  DATA, AND  NCSL  TRAJECTOR I ES; I F DESIRED 

APPLY  CNL  TO  INPUT 

READ( 5, 1 )MNVR 

WR I TE ( 6, 13 )MNVR 

DO  112  J«1,MNVR 

REAO(5 ,6 ,EN0”116)HPT , (NMA10( I, J), l “1,15) 

READ (5, 201 ) TITLE 

READ( 5, 688 9)  ( X ( K, J ) , K«1 ,NPT) 

IF( ICNL.NE.O)  CALL  NONLT(X(l, J) ,UNL,NPT) 

READ( 5, 201 ) TITLE 
READ(5,C889)(Y(K,J),K-1,NPT) 

CONTINUE 
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00560 
00570  113 
00580  C 
00590  C 
00600  C 
00610  C 
00620  C 
00630 
00640 
00650 
00660 
00670 
00680 
00690  C 
00700 
00710 
00720  105 
00730 
00740 
00750 
00760 
00770 
00780 
00790 
00800 
00810  111 
00820 
00830 
00840 
00850 
00860 
00870 
00880  115 
00890  116 
00900 
00910  1 
00920  2 
00930  3 
00940  4 
00950  5 
00960  6 
00970  7 
00980  8 
00990  9 
01000  11 
01010  12 
01020  13 
01030  15 
01040  48 
01050 
01060  201 
01070  202 
01080  203 
01090  204 
01100  415 
OHIO  416 


DO  113  K-1,NPT 
YPLOT(K,3)-0.0 

VERIFY  EACH  MODEL;  FOR  EACH  MANEUVER 


DO  116  M-l,NMOD 

WRITE(6,415) (NAMEMt I ,M), 1-1,15) 

WRI TE ( 6, 204 ) 

DO  115  J«1,MNVR 

WRITE(6,416)( NMA 10(1, J), 1-1, 15) 

CALL  RESP0N(YREC(1,J),X(1,J) , N,GAMMA1 ( 1 ,M) , XLAMDA, NPT) 

COMPARE  Y AND  YREC.  CAL  RMS  ERROR  ♦ *RMS  ERROR 
SUM-0.0 

DO  105  K-l.NPT 

SUM- SUM* YREC ( K,J)«YREC(K,J) 

RMS-SQRKSUM/NPT) 

SUM-0.0 
SUMW-0.0 
I »J 

DOlll  K-l, MPT 
YPLOT(K,l)-Y(K,l) 

YPLOT(K,2)=YREC(K, I ) 

SUMW-SUMW*Y(K, I )-Y(K, I ) 

SUM-SUM* (YU,  IO-YREC(K,l))*(Y(K,l)-YREC(K,l)) 

RMS-SQRKSUM/NPT) 

RMSW-SQRT ( SUMW/HPT) 

PC-100. -RMS/RMSW 
WRITE<6,48)  I, RMS, PC 

I F < I PLOT . NE . 0 ) CALL  PL0P8 (NPT, 3, YPLOT, 200 , . 2 , . 2,5,NAMEM( 1,M) , 

* NMA 10(1, J), BUFF) 

CONTINUE 

CONTINUE 

CALL  PICSIZtO. 0,0.0) 

FORMAT (4I5,F10.0) 

FORMAT ( 4F20 . 0 ) 

FORMAT (2I5,F10.0) 

FORMAT (/, 5X, ’ Z DOMAIN  DENOMINATOR') 

FORMAT ( /, 5X, ' Z DOMAIN  NUMERATOR') 

FORMAT( I 5, 15A4) 

FORMAT(//5X,15A4, 'POINTS  - ',15 ,//) 

FORMAK5X, 'THE  SYSTEM  IS', 13,'  BY  ',13,'  AND  ',13,'TH  ORDER') 
FORMAT ( 8F10 . 0 ) 

F0RMAT(//,5X, 'VERIFY  ',12,'  MODELS  OF  ORDER  ',12,//) 

FORMAT ( 10X, ' I CNL  - ' , I 2,/,10X, ' I PLOT  - ' , I 2, / , 10X, ’ I DY  - ',12) 
FORMAT ( 10X, 'MNVR  - ',12 ,///////) 

FORMAK15A4) 

FORHAT(//,7X,'RMS  ERROR  ',11,'  - ', 

• F10.5,5X,'%  RMS  ERROR  «',E1G.7,//) 

FORMAT ( 80A1) 

FORMAT (2X, ////////////////////) 

FORMAT (2X,80( )) 

FORMAT  C5X,40('*')) 

FORMAT(//,5X,15A4,//) 

FORMAT (//,2X,15A4,//) 
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TT 


m,  Pk0I  IS  BEST 

pbo,  co?x  yuaiisH®  io  oo.c  — ' 


01120  6089 

F0RMAT(5E16.7> 

01130  C 

01140 

STOP 

01150  C 

01160  C 

LIST  OF  VARIABLES 

01170  C 

01180  C 

N « ORDER  OF  MODEL 

01190  C 

NMOD  « NUMBER  OF  H(Z)  MODELS  (TO  BE  VERIFIED  WITH  SINGLE  CNL) 

01200  C 

E.G.,  10KN,15KN,&  20KN  MODELS  WITH  15  KN 

CNL 

01210  C 

NAMEM-  MODEL  TITLE  CARD 

01220  C 

NMA 1 0=  TITLE  OF  I/O  SET 

01230  C 

GAMMA=  PARAMETER  VECTOR  OF  A MODEL 

01240  C 

GAMMA 1*  STORAGE  ARRAY  IN  WHICH  ALL  MODELS  ARE  STORED 

01250  C 

1 CNL*  1 IF  CNL  USED,  0 OTHERWISE 

I 

01260  C 

1 PLOT*  1 IF  CALCOMP  PLOT  DESIRED,  0 OTHERWISE 

01270  C 

MNVR  « NUMBER  OF  INPUT  OUTPUT  SETS  TO  BE  VERIFIED 

01280  C 

X(K, J)*  INPUTS;  J*l,...,NIO 

01290  C 

Y(K,J)»  INPUTS;  J=1,...,MNVR 

01300  C 

YRECCK, J)*  NCSL  RESPONSES;  J-1,...MNVR 

01310  C 

UNL(K)*  NONLINEAR  TABLE 

01320  C 

RMS  « RMS  RECONSTRUCTION  ERROR 

01330  C 

PC  -PERCENT  TOTAL  ENERGY 

01340  C 

i * 

01350  C 

• 

t 

01360 

END 

r 

01370 

SUBROUTINE  PRVEC(A,N) 

• 

! 

01380  C 

THIS  SUBROUTINE  OUTPUTS  DOUBLE  PRECISION  SINGLE 

DIMENSIONED  ARRAY 

1 

01390 

DIMENSION  A( 1) 

j 

01400 

DOUBLE  PRECISION  A 

i 

01410 

WR 1 TE( 6, 31 ) 

01420 

WRITE (6,1)(A(I),I*1,N) 

01430  1 

FORMAT( IX, 5G20 . 10 ) 

1 . 

01440 

WR 1 TE ( 6, 31 ) 

t 

01450  31 

FORMAT (/) 

01460 

RETURN 

01470 

END  . 

i 

01480 

SUBROUTINE  NONLT(U,UNL,NPT) 

01490 

DIMENSION  U( 1 ) ,UNL( 1 ) 

01500 

0010  K*l, NPT 

01510 

UABS-A8S(U(IO) 

\ 

01520 

1 *UABS*1 .0 

j 

01530 

IP1-I+1 

01540 

UTEM-UNK 1 )*(UABS- 1+1 .0)*(UNL( 1 P1)-UNL( 1 ) ) 

01550 

U(K)-UTEM«U(K)/UABS 

01560  10 

CONTINUE 

I 

01570 

RETURN 

i 

01580 

END 

01590 

SUBROUTINE  RESPONC X, V, N, GAMMA, XLAMDA,MP1) 

01600 

DIMENSION  X(1),V(1), GAMMA ( 1 ) , XLAMDA ( 1 ) 

j 

01610 

REAL* 8 XSAV, GAMMA, XLAMDA 

I 

i 

* 

01620 

NM1-N-1 

i 

01630 

NP1-N+1 

01640 

NPNP1-N*N*1 

t 

01650 

NPNP2-N*N*2 

1 

01660 

00191*1, NPNP1 
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i 

i 

01670  19 

XLAMDAU >-0.0000  j 

01680 

XSAV-0 . ODOO  ' 

01690 

D020K«1,MP1  1 

01700 

IF(N.EQ.l)  G0T025 

* 

01710 

D021 1 -1, NM1  1 

• 

01720 

J-NP1-I 

01730  21 

XLAMDA( J) -XLAMDA( J-l) 

0171)0  25 

CONTINUE 

01750 

D022I«1,N 

01760 

J-NPNP2-I 

01770  22 

XLAMDA( J) =XLAMDA( J-l) 

01780 

XLAMDA( 1 ) =XSAV 

01790 

XLAMDA(NP1 ) »V( K) 

01800 

XSAV-0 .0000 

01810 

D023 1 “1 , NPNP1 

01820  23 

XSAV-XSAV- GAMMA ( 1 + 1 ) *XLAMDA( 1 ) 

01830 

1 F( DABS (XSAV) .GE.1.0D10 ) XSAV-0. ODOO 

01840  20 

X (K) -XSAV 

01850 

RETURN  . - | 

01860 

END 

01870 
01880  C 

SUBROUTINE  DASH(X,Y,N, INC,DL,SL) 

• 

01890  C 

THIS  SUBROUT  1 HE  DRAWS  A DASHED  LINE  DEFINED  BY  THE  ARRAYS  X AND  Y. 

01900  C 

X IS  THE  ARRAY! 1)  CONTAINING  THE  X COORDINATES. 

01910  C 

Y IS  THE  ARRAY ( 1 ) CONTAINING  THE  Y COORDINATES. 

01920  C 

N IS  THE  HUMBER  OF  POINTS  IN  X AND  Y. 

i 

01930  C 

INC  IS  THE  INCREMENT  (CAN  BE  NEGATIVE,  X AND  Y WILL  THEN  USE  NEGATI 

01940  C 

ARGUMENTS,  BUT  N IS  STILL  POSITIVE). 

01950  C 

DL  IS  THE  LENGTH  OF  THE  DASH. 

J 

01960  C 
01970  C 

SL  IS  THE  LENGTH  OF  THE  SPACE  BETWEEN  DASHES. 

* 

01980 

REAL  X(1),Y(1),R(2),P(20),Q(20) 

01990 

IF(N)6,5,4 

02000  4 

IF(N-1)5,5,7 

s 

07010  5 

CALL  PLOT ( X, Y,N) 

• 

02020  6 

RETURN 

02030  7 

L-l 

02040 

K-IABS(INC) 

i 

02050 

M-l 

02060 

A-0.0 

02070 

R( 1 ) -DL 

02080 

R( 2 ) -SL 

02090 

PCl)-X(l) 

i 

02100 

a<l>-Y(l) 

02110 

1-1 

02120 

DO  38  J-2,N,K 

02130 

L-L+INC 

1 

02140  18 

1-1*1 

‘ 

02150 

P(  1 )-X(L) 

02160 

Q( 1 ) *Y( L > 

02170 

B-SQRTC  ( P(  1 )-P(l-l))««2*(Q(l )-Q{ 1 -1 > )**2)*A 

02180 

IF(B-RCM) >33,26,23 

i 

02190  23 

A»(R(M)-A)/(B-A) 

j 

02200 

P( 1 )»(P( 1 >-P( l-l))*A*P( 1-1) 

02210 

Q( 1 )-(Q( 1 )-Q( 1-1) )*A*Q(  1 -1) 

02220  26 

A-O.ODO 
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02230 

1 F (M . EQ. 1)CALL  PLOT(P,Q,l> 

0221*0 

M-3-M 

» 

02230 

P(l)-P(l) 

02260 

QCl)-Q(l) 

02270 

I«1 

02280 

GO  TO  18 

02290 

33 

IF(l-20)38, 34,34 

02300 

34 

IF(M.EQ.1)CALL  PLOT(P,Q,l) 

02310 

P(l>-P(l) 

02320 

Q(l)-Q(l) 

02330 

1 -1 

02340 

38 

A«B 

„ 

02350 

IFC 1 . NE.l. AND, M.EQ.l) CALL  PLOT(P,Q, 1) 

02360 

GO  TO  6 

02370 

END 

02380 

SUBROUTINE  MYAX 1 S( XO, XF , YO, YF,NLD, NSD,T 1 C) 

02390 

C 

02400 

C 

THIS  SUBROUTINE  PLOTS  A TICKED  LINE  ON  THE  XY  PLOTTER  FROM  (XO,YO) 

02410 

C 

(XF,YF)  WITH  NLO  LARGE  DIVISIONS  AND  NSD  SMALL  DIVISIONS  PER  LARGE 

02420 

C 

WITH  A TIC  OF  SIZE  TIC. 

02430 

C 

02440 

REAL  X(3),Y.S) 

02450 

OX"(XF-XO)/(NLD*NSD) 

02460 

OY*(YF-YO)/(NLD*NSD) 

. , 

02470 

A*SQRT(DX**2+DY**2) 

02480 

DXT— TIC*DY/A  ’ 

. 1 

02490 

DYT«TIC*DX/A 

1 

02500 

X(3)*XO 

1 

02510 

Y(3)-YO 

i 

02520 

DO  20  1-1, NLD 

"« I 

02530 

B-1.0 

02540 

DO  20  J"1 , NSD 

.! 

02550 

X ( 1) *X ( 3 ) +B*DXT 

! 

02560 

X(2)-X(3) 

• 

02570 

X(J)mX(3)*DX 

1 

02580 

Y ( 1) «Y ( 3)+B*DYT 

.1 

. j 

02590 

Y(2)-Y(3) 

• 

02600 

Y( 3 )“YC  3 ) +DY 

I 

02610 

CALL  PLOT(X,Y,3) 

• | 

02620 

20 

B-0.5 

02630 

X ( 2 ) «X( 3 ) + DXT 

i 

02640 

Y( 2) *Y( 3 )+DYT 

i 

02650 

CALL  PLOT (X(3),Y(3),2,-1) 

t 

02660 

RETURN 

02670 

• END 

02680 

SUBROUTINE  SCALE(ND,NP,NF,X,NX,XO,XF) 

02690 

C 

» 

02700 

C 

THIS  SUBROUTINE  FINDS  QUANTIZED  UPPER  AND  LOWER  LIMITS  OF  A DATA  SET, 

! 

02710 

C 

INPUT: 

.# 

02720 

C 

IN  CHOOSING  SCALES  FOR  GRAPHS. 

. ♦ 

•i 

02730 

C 

ND  IS  THE  NUMBER  OF  DIVISIONS  (QUANTUM  LEVELS). 

i 

! 

02740 

C 

NP  IS  THE  NUMBER  OF  POINTS  PER  FUNCTION. 

. \ 

02750 

C 

NF  IS  THE  NUMBER  OF  FUNCTIONS. 

1 

02760 

C 

X IS  THE  ARRAY(NX.NF)  CONTAINING  THE  FUNCTIONS. 

♦ 

♦ • 

02770 

C 

NX  IS  THE  DIMENSION  OF  X. 

. I 

02780 

C 

OUTPUT: 

•i 

t 


02790  C 
02800  C 
02810  C 
02820 
02830 
02840 
02850 
02860 
02870 
02880  8 
02890 
02900 
02910 
02920 
02930  13 
02940 
02950  15 
02960 
02970 
02980 
02990 
03000  20 
03010  21 
03020 
03030 
03040 
03050 
030C0 
03070 
03080  C 


.1 

03090  C 

THIS  SUBROUTINE  PLOTS  ONE  OR  MORE  FUNCTIONS  OF  A CONSTANT  INTERVAL  ON 

vi 

03100  C 

X-Y  PLOTTER. 

;•! 

••  1 

03110  C 

NP 

IS  THE  NUMBER  OF  POINTS  PER  FUNCTION. 

i 

1 

03120  C 

NF 

IS  THE  NUMBER  OF  FUNCTIONS. 

• i ' 

03130  C 

DATA 

IS  THE  ARRAY(ND,HF ) CONTAINING  THE  FUNCTIONS. 

1 

03140  C 

ND 

IS  THE  DIMENSION  OF  DATA. 

. 1 

. 

03150  C 

TO 

IS  THE  LOWEST  VALUE  OF  THE  INDEPENDENT  VARIABLE. 

03160  C 

T 

IS  THE  INTERVAL  BETWEEN  TWO  SAMPLE  POINTS. 

I 

03170  C 

NSZ 

IS  THE  HUMBER  OF  PLOTS  PER  PAGE  (SIZE  CONTROL). 

03180  C 

LABEL 

IS  A HOLLERITH  FIELD  TO  BE  WRITTEN  ACROSS  THE  TOP  OF  THE 

03190  C 

(MUST  END  WITH  A ) , 

■ \ 

03200  C 

INDEEP 

IS  A HOLLERITH  FEILD  DESCRIBING  THE  INDEPENDENT  VARIABLE 

i 

03210  C 

END  WITH  A ). 

\ 

03220  C 

BUFF 

IS  A WORKING  ARRAY(NP,NF+1) 

03230  C 

• i 

03240 

REAL  DATA( ND, NF ),TO,T,BUFF(NP,1),SLEFT/0.0/ 

j 

03250 

INTEGER  LABEL(l), IN0EEP(l),NSIZE/0/ 

| 

03260 

S-1.0/NSZ 

.03270 

IF(SLEFT. 

GE.S.ANO.NSZ.GE.NSIZE)GO  TO  9 

03280 

CALL  PICSIZ(16.0*S,28.0) 

03290 

SLEFT-1.0 

* 

03300 

NSIZE-NSZ 

03310  9 

CALL  ORIGIH(14.0»S,28.0*(1.0-SLEFT)*S*2.5) 

i 

03320 

CALL  FACTORC1. 0/S, 1.0/S) 

■ j | 

03330 

CALL  MYAXIS(-12.0,0.0,0.0,0.0,12,5,-.2S) 

03340 

CALL  MY AX  ISC  0.0, 0.0, 0.0, 20. 0,10,5,-. 5) 

THIS  PAGE  IS  BEST  QUALITY  FRACTICABL* 
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XO  IS  THE  LOWER  LIMIT. 

XF  IS  THE  UPPER  LIMIT. 

REAL  X(NX,NF),XS(U)/1.,1.2S,1.5,2.,2.5,3.,4.,5.,6.,8.,10./ 
A*X( 1, 1 ) 

B«A 

00  8 l-l,NF 
DO  8 J-1,NP 

IF(A.GT.X(J,I))A*X(J,I) 

IF(B.LT.X(J,I))B-X(J,I) 

C=0.1*(ND-1) 

B-B-A 

I F < B . EQ. 0 . 0 ) B = C 
I F ( B . LE . C ) GO  TO  15 
C-10.0*C 

I F ( B . GT . C)GO  TO  13 
C»0.1*C 

I F ( B . LE . C)GO  TO  15 
B = B/C 

00  20  1-1,11 
IF(B.LE.XS( I ))G0  TO  21 
CONTINUE 
B«XS(I)*C/(MD-1) 

X0“B* I F I X(A/B) 

IF(A.LT.O.O)XO»XO-B 

XF*XO+ND*B 

RETURN  • 

END 

SUBROUTINE  PLOP8 (NP,NF, DATA, ND, TO, T,NSZ, LABEL, INDEEP, BUFF) 
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03350 

CALL  MYAX ISC  0.0,-12.0,20. 0,20. 0,12,5,-. 25) 

05560 

CALL  MYAXIS(-12. 0,-12. 0,20.0, 0.0,10, 5, -.5) 

03370 

CALL  SCALE(12,NP,NF,DATA,ND,A,B) 

05380 

C=( B-A) / 12 . 0 

03390 

DO  19  1 =1, NF 

03400 

DO  19  J=1,NP 

03410 

19 

BUFF ( J, 1 )=(A-DATA(J, 1 ))/C 

03420 

N=NF+ 1 

05430 

A=20 . 0/ (NP-1 ) 

03440 

DO  23  1=1, NP 

03450 

23 

BUFFC 1 ,N)=( l-l)*A 

03460 

DO  26  1=1,13 

03470 

CALL  NUMBERC 1-12.825,- 2.375, .25*S,B, 90. 0,3) 

03480 

B = B-C 

03490 

26 

IFCABS(B) .LE.1.0E-5*ABS(C) )B=0.C 

03500 

A= ( NP-1 ) *T*  0 . 1 

03510 

B=T0- A 

03520 

DO  29  1=1,11 

03530 

29 

CALL  NUMBER ( .85,2.0*I-2.2,.25*S,I*A+B,90.0,3) 

03540 

CALL  SYMBOL ( 1. 6, 4. 0, ,5*S, INDEEP, 90.0) 

03550 

CALL  SYMBOL (-12. 75, -2. 375, .4*S, LABEL, 90.0) 

03560 

CALL  PLOT( BUFF (NP, 1 ) , BUFF(NP,N) ,NP, -1) 

03570 

1 F(NF . EQ. 1 )GO  TO  46 

05580 

A = 0 , 1 

03590 

B = 0 . 2 

03500 

1 = 1 

03610 

J = 1 

03620 

DO  45  K=2,NF 

03630 

CALL  DASH ( BUFF (U,K), BUFF ( J, N), NP, l,B, A) 

03640 

8 = 0 . 5 *B 

C3650 

1 F ( B . GE . A) GO  TO  44 

03660 

A = 2 . 0 *A 

0367C 

B = 2 . 0 *A 

03680 

44 

l=-l 

036S0 

45 

J=NP- J+l 

03700 

46 

SLEFT=SLEFT-S 

03710 

CALL  FACT0R(1. 0,1.0) 

03720 

RETURN 

03730 

END 

03740 

//GO. 

PLOTMSG  DD  SYSOUT=A 

03750 

//GO. 

PLOTTER  DO  SYSOUT=P 

03760 

//GO. 

SYS  IN  DD  » 

03770 

2 

1 

03780 

1 

1 

037S0 

5 

03800 

// 

END  OF 

: DATA 
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